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AN IMPROVED METHOD OF NEUTRON SPECTROSCOPY 
USING THRESHOLD DETECTORS 


Arthur Degen Kohler, Ir. 


Lawrence Radiation Laboratory 
University of California 
Berkeley, California 


December 11, 1964 


ABSTRACT 


An improved technique of measuring neutron~energy spectra 
in the range 2.5 to 30 MeV has been developed. In this technique 
threshold detectors are exposed to the unknown neutron flux, The ac- 
tivities of residual nuclei from (n,p), (n,a), and (n, 2n) threshold re- 
actions are measured by gamma-spectrum analysis. Experimental 
excitation functions are used, Measured activities and cross sections 
form a set of integral equations which are solved for the neutron 
spectrum by using a least-squares technique, 

Neutron-reaction cross sections are not well known above 
about 20 MeV, and 30 MeV is selected as the highest energy for 
which this method is feasible based on available data. An approximate 
correction is made to measured activities to account for activation 
caused by neutrons with energies above 30 MeV. 

The previous technique, on which this work is based, was im- 


proved by: 


(a) substitution of experimental reaction cross sections (from re- 
cently published literature) for those calculated from the continuum 


model of the nucleus; 


* 
Lieutenant, Civil Engineer Corps, United States Navy. 














(b) solution of the spectrum from the set of integral equations by 


a least-squares technique; 
(c) calculation of the amount of activation of threshold detectors 


caused by neutrons with energies above 30 MeV within a thick concrete 


shield exposed to a 6.2-BeV proton beam, 


Application of this improved technique is demonstrated in ex- 
periments with a Be (d, n)B*° neutron source, and the neutron spec- 


tra in a thick concrete shield exposed to a 6.2-BeV proton beam, 
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I, INTRODUCTION 


Activation of threshold detectors by (n, p), (n,a) and (n, 2n) 
reactions has been developed and. demonstrated as a useful method of 
neutron spectroscopy ina previous study. : This method is applicable 
to the energy range of 2 to 30 MeV andis sufficiently sensitive to 
measure very small fluxes of the order of 1 to 100 n/cm* sec, 

In this method, a number of different detectors, usually in the 
form of 4-in. diam metal discs are exposed to the unknown neutron 
flux. After exposure, the detectors are moved to a low-background 
counting facility. The gamma-ray spectrum of the activated detector 
is measured by using a Nal crystal, photomultiplier, and pulse-height 
analyser. Analysis of the gamma-ray spectrum identifies the radio- 
active nuclei and their decay rate. In application, a set of metal-disc 
detectors of aluminum, iron, cobalt, and nickel might be exposed to 


the neutron flux. In the aluminum disc one might observe gamma rays 
27 27 27 
of Mg 


and Nace which were created in Al" (n, p)Mg- and Al 
(n, eiNac* reactions, The half lives of residual nuclei must be suffi- 
ciently long, in our technique, to allow the detector to be moved from 
the neutron flux to the gamma-ray counter. Our gamma-ray measure- 
ments have benefited from improvements in the low-background count- 
ing facility. A 
The relation between detector activity, cross section, and 


neutron flux may be expressed as 


foe] 
0 
where A' = detector saturation activity or disintegration rate 


‘ -41 
per target nucleus in sec |, 


2 -1 4 


o(E)= neutron flux incm “sec “MeV, 





: . ; 2 
o{E) = reaction cross sectionincm , 


E = neutron energy in MeV, F 
i = subscript identifying the reaction and residual 
nucleus, 


There is one such equation for each reaction observed, A number of 
reactions form a set of integral equations which, under favorable con- 
ditions, may be solved for the unknown flux. Our previous solution 
method was a matrix inversion in which the form of 6(E) was assumed, 
and numerical values of 6(E) were found. This solution, when sub- 
stituted back into the equation, reproduced exactly the activities, A.- 
In cases where the original activity measurements contain substantial 
experimental errors, the flux solution may become unstable and pro- 
duce negative values of o(E), which are meaningless, We encountered 
this difficulty and substituted a solution method which solves by iter- 
ation, using only positive trial fluxes, and seeks to produce a set of 
calculated activities, based on the trial flux, which best fit the meas- 
ured set of activities. 


It was desirable to incorporate recent experimental threshold 


alae ribs wt a ted ap aah a in) niga a vag con aa! see uae Sy a da Rp ER ae cae ast 


reaction cross sections (excitation functions) inthe method, Our pre- | 
vious calculations of reaction cross sections did not show satisfactory 
agreement with recently published experimental data, The calculated 


cross sections were discarded in favor of experimental data, 
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II. MEASUREMENT OF DISINTEGRATION RATE 


A. Counting Equipment 


The counting equipment is similar to that in our earlier report 
on threshold detectors. : A single 4-in. -diam by 2-in. -thick Nal(T1 ) 
scintillation crystal is coupled to a DuMont type 6363 photomultiplier 
{PM) tube. The crystal and-PM tube are located in a new large, low- 
activity concrete cave, : Inside the cave a smaller 2-in. -thick lead 
shield surrounds the crystal and PM tube. Within the y-ray energy 
snterval 0,080 to 2.080 MeV the background count rate is 141.6 counts/ 
min and is quite constant, The previous setup used a second crystal 
‘or continuous background measurement, and this crystal is no longer 
necessary. | . 

A Penco PA-4 100-channel differential pulse-height analyzer 
(PHA) is used to analyze and store the y-ray data. Gain is set to 
cover the energy interval 0,080 to 2.080 MeV in 100 channels, each 
20-keV wide. Stored PHA data is tabulated on a Victor printer adding 
machine and the counts/channel vs channel number are simultaneously 


plotted on a Mosely Autograf, 


B. Gamma-Ray Spectrum Analysis 

1. Photopeak Analysis 

Our method of y-ray spectrum analysis is essentially the same 
2s previously reported and is briefly outlined here. , Our method isa 
variation of spectrum stripping. Figure 1 shows one of our measured 
ycray spectra from an Al detector. Data points are shown as dots. 
The prominent peak near channel 60 identifies Na?*, The smaller 
peak near channel 86 represents a sum peak for scintillations from the. 
annihilation y ray andthe 1,28-~MeV y ray in coincidence, It is of no 
znterest in this analysis. Figure 1 also shows this spectrum separated 


tnto its components, We are interested in calculating the photopeak 


area in counts/sec, Analysis proceeds as follows: 
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Fig, 4. Gamma-ray spectrum including Na“ and 
its components, 
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(a) Background is subtracted from the spectrum, 

(>) The response of unidentified gamma rays of higher energy than 
t2e photopeak of interest elevates the peak. A constant quantity, Y, 
a= estimate of this response, is subtracted from the spectrum, 

(c) Based on the position and height of the photopeak maximum the 
Compton distribution is found from a library of such functions and 
ecitracted from the spectrum. 

(d) What remains, in principle, is the photopeak only. Aside from 
statistical error or errors in the quantities subtracted (a,b, and c 
above), the photopeak is Gaussian in shape. This remaining data is 
Stted to a Gaussian curve by a nonlinear least-squares technique which 
selects best values of photopeak area, variance, and mean (position of 
peak maximum), 

(e) Photopeak area is the number of counts in the peak per time. 
Cerrections are made for PHA block time, y-ray self-absorption in 
the detector, branching ratio for the y ray observed, and the ratio of 
counts in the peak to y rays emitted by the activated detector, The 
result is the disintegration rate of the isotope of interest in the detec- 
sor. The calculation is coded in FORTRAN for the IBM 7090 computer 


+= &@ main program titled SUPER-3 and 15 subroutines. 


2. Least-Squares Analysis 

Our method of y-ray analysis uses data only from the photo- 
peak, which is 5 to 15 channels wide out of the 100 channels of data. 
A newer method currently in great favor is least-squares analysis,’ 
wich uses data from most of the spectrum, This method assembles 
a composite spectrum which is the sum of individual spectra of gamma 
rays of known radioisotopes. The heights (which are proportional to 
t£e activities) of the individual spectra are adjusted to minimize the 
sum of the squares of the differences, channel by channel, between the 


composite spectrum and the experimental spectrum of unknown com- 


position. 











Least-squares analysis should be a more accurate method thar. 


photopeak analysis. There are several practical reasons why this has 
not been attempted here. Many of the y-ray spectra we have analyzeé 
to date were very complex, being the result of exposing detectors to 
particle radiation whose energy spectrum reached very high energies 
(for example, the 6.2-BeV Bevatron), High-energy particle radiation 
can produce very many different radionuclides. Least-squares anal- 
ysis requires a library of y-ray spectra, representing each radionu- 
clide likely to be encountered. Further the spectrum of a single 
radionuclide is seen to be different with different detector materials 
or different detector thicknesses, (All our detectors are "thick" 
sources, between 1/32 and iin.). For our use to date, applying 
least-squares analysis would be a monumental task. Where thresh- 
old detectors are standardized to several materials, each of a stand- 
ard thickness, and exposed to particle radiation not having a signifi- 


cant ultra-high-energy component, application of least-squares gam 


analysis should be practical and advantageous. 
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Il. THRESHOLD DETECTORS 


A. General Requirements 


After some experience, the following criteria guide our 


#e.ection of threshold detectors: 


a. The Q value* for the reaction must be negative or, if not, 
t2ere must be an-effective energy threshold, so that neutrons below 
sal energy do not initiate the reaction. 

=. The material must be available in sufficient purity that activa- 
t:om of impurities does not obscure the reaction of interest. 

<. The material must be safe to use and reasonably inexpensive. 

¢. The residual nucleus must be radioactive and emit a y ray with 
4+ :centifiable photopeak. 

e. The half life of the residual nucleus must be of a convenient 
-ezgth. The shortest half life we have positively observed is 9.5 min 


Ls mMe?’, The longest is 5.27y for er 


after about one year follow- 
-=g :rradiation, It is preferable to have half lives between 1h and 70 
4a7s. 

‘. The reaction cross section must be an experimentally known 
feaction of neutron energy. This last criterion largely limits our 
s&a:ce to (n, p), (n,a), and . (n, 2n) reactions, Other experimenters 


5,6 


save also used (n,n!') reactions. 


Our previous study led us to calculate cross sections by using 


, d 4 
s2£{e continuum model of the nucleus. s158 


This was largely because of 
a= absence of sufficient experimental data. Recent experimental data 
save Degun to fill in energy gaps and improve resolution of a number 


: reactions, Figures 2 and 3 show examples of relatively good and 


al 





The Qvalue for the reaction X (a, b)Y is defined as the energy equivalent 


3: the mass difference between the input components, X plus a, and 


t2e reaction products, b plus Y. 
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Fig. 3. Cross section vs energy for Go? fa, a)Mn-~, 
The dashed line shows the result of the continuum- 
theory calculation, 
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poor agreement, respectively, between experimental measurements 
and our calculated cross sections. Experimental data is indicated by 
dots, and no attempt is made to identify experimenters, as the data 


9-33 


comes from many sources, More refined continuum-theory cal- 


culations show better agreement but are still considered inadequate 


: 20,241, 24,25,28 
for our requirements. 


B. Useful Reactions and Cross Sections 


Table I lists reactions which were found useful in our exper- 
iments, and some of their properties. Table Ii lists a number of 
other reactions of interest in future experiments, for which exper- 
imental cross sections are determined or are being determined, 
These lists do not exhaust all possibilities but contain a variety of re- 
actions which should satisfy many requirements. Figures 4 through 
13 show experimental-cross-section measurements (dots) taken from 
the literature and excitation curves (solid lines) derived from these 
data, The excitation function is a continuous function describing the 
variation of the reaction-cross section with energy. We have approx- 
imated excitation functions by smooth-drawn curves through the exper- 
imental points, Increased weight was given to newer data. Very 
little experimental data exists beyond 21 MeV. Between the last data 
points and 30 MeV, the excitation curves are extrapolated. The shape 
drawn in this region.is partly inspired by proton cross sections at 
high energies (see Figs, 23-27 of this report), 

All detectors used in our experiments were in naturally occur- 
ring isotopic abundances, As we intend to solve neutron spectra up 
to 30 MeV and to correct for activation by higher-energy neutrons, we 
must account for the several isotopes producing the same residual 
nucleus through different reactions, For example, in natural nickel, 


58. . ‘ 
neutrons. may produce Co” in a variety of reactions: 
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+He8t pic abundance Reaction Threshold energy 
(73 pes (MeV) 
67.8 Ni tapiGe” 0.4 
26.2 NiO? (n, t)C0°? 14.7 
i 
Similarly neutrons en ie and nio4 can produce’Co”®: i 


£a:.tat:on functions represented in this report include reactions in- 
@a¢<2 in all natural isotopes. Experimental measurements of cross 
ee;ti92s are often made with natural elements, so this is not a serious 


oes em. 
8 


AS an example, the excitation curve for Fig. 4, labeled Ni? 
5 58 


tt. o*C_”, is actually intended to represent the cross section for Co 


ges action in natural nickel per Nin? nucleus. 
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Table I. Threshold detectors. 





Reaction ‘Theoretical Effective Half life of Form of — 
threshold threshold residual detector 
(MeV) (MeV) nucleus 
Ni??(n, p)Co?? -0.4 4.2 71 days 4-in, metal disc fie 
Ain pio" 1.8 Oe ‘956imnin. -4inGHelal disc 
| Co? {n, ayMn?° -0.3 5.3 2.58 h. — 4-in, metal disc ME 
| Fe? °(n, p)Mn2° 2.9 5.0 2.58 h 4-in, metal dise 
Ti*8(n, p)Se*? 3.2 5.2 44,0 h 4 anv tetal, diese 
Mg~*(n, p)NaZ* 4.7 6.4 15.0 h 4-in, metal dise 
Al" (n, a)Na“* 3.4 5.9 45.0 h Ain. metal dise 
ace ry ea 9.3 9.4 13.2 days Boxed crystals 
Co”’(n, 2n)Co>® 10.2 40.8 71 days © 4-in. metal dise! 
Ni? Fin, anya?’ 14.8 12.5 36h 4-in, metal disci 
“The theoretical threshold is calculated as -Q X de where Q is 


the Q value for the reaction and M is the mass number of the targe s 
-nucleus. The effective threshold is the energy at which the cross 


section is 1/100 of.its peak value. 
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Table II, . Additional threshold reactions 
eS eee 





Reactions Theoretical Half life of residual : 
threshold nucleus : 
. (MeV) } 
ea * =. p)Gu"* 0,2 42.9 hr : 
a i 
25° in, p)Fe?? 0.8 45 days | 
38s 52 
+2 “(s, p)V 3.1 3.77 min | 
Ms (2,a)V°7 0.6 3.77 min 
ag 48 
Hd “is, a)Se 1.6 44,0h 
were 202 
ae {n, 2n)T1 8.8 42.0 days 
q*¢ 
:°- "tn, 2njsp'2° 8 9.3 16 min 
ed OME Tor can 10-6. 6S ais | 
<>": Pajeul* 40.4 12.9 h : 
2... 5 44m i 
fe “in, 2n)Sc 11.6 2.4 days : 
kz "(x, 2n)Na“” 42.9 2.6 y : 
“<=- i 
: F°*!e, onypt8 44.0 1.87 h ‘ 
Dy is : te"4 5, 2nyzn®? 41.8 38 min | 
target re n, 2n)Fe-? 13.6 9.0 min 1 
s r’ "is, 2nyp?? 12.7 2.55 min 1 
aq i 
9°" 7(n, balin’ @ ec 9.4 50.0 days 
Aa ’(n, 2n)As/4 10.4 18 days : 
¢37 196 g . 
wa" (nn, 2n)Au 8.4 6.1 days : 
id Ue ete 20.3 20.4 min ' 
a 28 ; 
& "tn, p)Al 3.99 2.3 min 
mim | 
$< taapini 4.3 2.65 h | 
ere : 
2? * Fn ntyntt> ™ 0.2 4.6 h 
as*? Fen, n'yRn is m 0,04 54 min 
74 
=e eeyaie® 2.0 a aan 


a? 
Ms “(n, 2n)MoX\ 13.3 46 min 
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Fig. 6, Excitation function for Co” (n, a)Mn° 
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Fig. 9, Excitation function for Mg~*(n, p)Na“*, 
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IV, NEUTRON-ENERGY SPEC TRUM 


A. Basic Relationships 
eps. 


The measured detector activity may be reduced to saturation § 
activity (see Appendix A) and this, in turn, related to cross section 


and neutron flux by 


ioe) 


A! = f $(E) o,(E) dk, (2) 
0 
where A,' = detector saturation activity or disintegration rate per g 


target nucleus in sec , 


o(E)= neutron flux in ea sec”? Mev"? 

o(E) = reaction cross section in arg 

E = neutron energy in MeV, and 

i = subscript identifying the reaction and residual nucleuse 


There is one.such equation for each reaction observed. We wish to 
determine the neutron spectrum below 30 MeV, so the set of equations 


(2) is rewritten 


30 20 
A! = { $(E) o,(E) dE + | o(E) 0; (E) dE (3) 
0 30 
and 
A,! = A30, + Amax, : (4) 


Here we assume Dh gp is a relatively small part of A! and can be 
estimated separately (see paragraph V,B.1). Then we wish to solve 
the following set for o(E): 


30 
Asy = ‘i (E) o,(E) dE, (5) 
1 
0 





LoS 


Ev.atson (5) is a linear integral equation, called a Fredholm equation 


*2 t2¢ Inrst kind. General discussion of its solution may be found 


v-geuhere, 35 The following is a discussion of several practical solu- 
basen Sam smnethods to our particular problem, 
>ction 
B. Average Flux 
As an introduction, assume that we will solve for the average 
(2) +s& inthe interval 0 to 30 MeV. Any one of the set of equations (5) 
-¥ #sii:c:ent for this purpose, Rearranging 
rate per 30 a 
Aso : J o(E) o (E) dE =o ,x o_(E) dE, (6) 
0 0 
wer Pave 
ra: 
nucleu Ox, oe 30 : (7) 
sh to J a (E) dE 
quations 0 
where on is the average flux in cm74 Bee Mev !. However, the 
talse of On calculated will depend on the reaction chosen, as illus- 
(3) oF trated in Fig. 14. 
Figure 15 illustrates that the average -flux calculation may be 
*«',;ne< by defining an energy range over which the cross section is 
(4) “set ellective. The average flux calculated by Eq. (7) is limited to 
- "2 yange over which that cross section is effective. Using a number 
can be §£ +4 Selectors in this manner, one might make a crude estimate of a 


solve ; teectrum. A slightly more elegant method is considered next, 
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Fig. 15. Refined average flux calculation, 
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C. Step-Function Approximation 


In this method the flux spectrum is assumed to have the form 
of a number of level steps. See Fig. 16. The solution of this method 


is the calculation of the level of each step. 


4. Matrix Inversion 

For the case where the number.of flux steps, n, is equal to the 
number of reactions, and each reaction is unique or has a unique crosg 
section shape, the basic equation (5) relating activity, flux, and cross 


section may be written 
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To simplify notation, we define 


j 
I = i o_ (E)dE, (9) 


In matrix notation, Eq. (6) may be written 


[Azo] = [lo] IH, (40) 


and the. solution is 


[Azo]. 

This problem is simply to solve n linear algebraic equations in n 
unknowns, It is straightforward except that care is required in choosiy 
the energy array or location of energy steps between 0 and 30 MeV, 
These steps should be chosen to maximize the absolute value of determi 
nant [I], in order to insure accuracy in its inversion, [I] 4A full 
discussion of this feature is found in Ringle': report, 

A computer program has been written to select the best energy 
array from a number of possibilities and to solve for the step-function 
flux. Tests of this program have been constructed in which the ''true” q 
flux is chosen and the activities calculated from equation set (5), Such ¥ 
a test of the program is to solve for the flux steps, using these calcu- 
lated activities and the known cross sections, Our tests demonstrate a| 
program that satisfactorily reproduces the trial fluxes, Examples of 
such tests may be seen in Ringle's report. 

Our application of this method to experiments demonstrates a 
practical difficulty. In almost every case, one or more of the flux 
steps is negative and comparable in magnitude to the positive steps. 

The explanation is that the activities, Aj; are measured quantities andy 
subject to experimental error. If the true flux and cross sections were] 


known exactly, calculated activities would not exactly match measured ; : 
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#<t:vities because of experimental error, In order to reproduce 
s-easured activities exactly, the calculated flux often takes on negative 
taisces. This difficulty has been of such magnitude that matrix inver- 


w2on has been abandoned as a practical method, 


=. Least-Squares Fitting 

In this method, the number of flux steps, n, is chosen and 
ted not equal the number of reactions, £, 

The basic relation between activity, flux, and cross section is 
#ga:n Eq. (8). Suppose that a set of trial values are selected for the 


t-2x steps, Calculated "activities'! may be found from 


n 
B -) tri, (414) 
m ty 
j=4 
where Ba = calculated activity for reaction m, 
. = trial flux step between er and Es . 


vie degree of misfit between measured activities, A, }, and cal- 
m 


smeated activities, Bea}, may be computed as 


2 
y (so, 7 Pm] 
D = . (12) 
A 
m=1 


2 
30 
m 


The better the set of trial fluxes, if; + is chosen, the better the 
+a.c ulated activities, UB }, match the measured activities, Ao i, 
acd the smaller D becomes, An unstable solution is to be avoided by 
‘=o9 Sing only positive values for trial fluxes, The problem is to min- 


“=:ze D with respect to n variables, tf}. 


D. Polygonal Approximation 
teen Pproximation 


To obtain a flux-spectrum solution that is a continuous function, 


“sm cay assume that the solution has a polygonal form, composed of 


#*ta:ght-line segments connected at the ends (see Fig.17), The 
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Polygonal approximation of flux, 


Fig. 17, 
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solution of this method is the calculation of the flux at each line-seg- 
ment end point. Where the number of variables, 03 >4> P> cee , iS 
equal to the number of reactions with unique cross-section shapes, the 
spectrum may be solved by matrix inversion, - The procedure is 
similar to the matrix inversion previously applied to the step function 
and 1s subject to the same practical difficulty. The calculated flux os- 
tillates between positive and negative values when the activation data 
contain experimental errors. Least-squares fitting of positive values 
& ain provides the solution. 


The polygonal flux is assumed to have the form: 


E-E, E,-E 
(EF) = 65(e—--] + 9, (EOE (13) 
2 E, E, 4 E, Ey 
‘or E,< E<E,. Similarly, flux may be expressed between any other 


4 2 
two end points. The basic equation (5) relating activity, flux, and 


cross section may be written 





0 0 4 
ED E-E, \ E,-E 
+f >> E -E, + by EE oC (E)dE + 
/ 2 4 
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E-E 
n-4 30-E 
+f °n \30-E ony 3S :} 0 (E) dE, 
E n-1 n-1 


where n+ 4 is the number of variables, This equation may be further 


simplified to 
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Notation has been simplified by defining 
* 
‘oe =| Eo _(E)dE (45) 
mj m 
Ee 


and by using Eq. (9) to define Lay 


The number of variables, n+ 1, is chosen and may be less thar. 
or equal to the number of reactions, £. A set of trial values, {f. }, 
may be selected to represent the set of flux variables, 1; }. Based ox 


the trial values, calculated ''activities''’ may be found from 


n n-1 Z 


Sule e. ba : 
Ba) & Pow nj Bet ti) >, fhe. (Ea wereren | 
jt J j-1 jc ee 

(16) 
where Bun is the calculated activity for reaction m, Using Eq, (12) 
compute and minimize the degree of misfit between measured activities 
Aso }, and calculated activities, Bot by choosing better and better 
trial values, tf} Negative flux values are avoided by choosing only 


positive values for the trial variables. 
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E. Computer Programs 


Two computer programs were written to calculate flux spec- 
tra--FLUXSPOS for the step function and .FLUXPPOS for the poly- 
gonal flux. Least-squares fitting of trial fluxes to minimize D [see 
Eq. (12)] is carried out by subroutine VARMINT. ee These programs 
also calculate detector saturation activity per target nucleus caused by 
neutrons of less than 30 MeV; refer to Eq. (3) and the Appendix, Input 
data for the programs includes disintegration rate (an output of SUPER-3), 
count time, wait time, irradiation time, half life, detector weight, 
atomic weight of target element, abundance of target isotope, fraction 
of activation caused by neutrons below 30 MeV, neutron-energy end 


points, and initial trial values of flux, Calculations are done on the 


Lawrence Radiation Laboratory's IBM 7094 computer. 























V. EXPERIMENTAL RESULTS 


A. Measurement of Be (d, n) Bo” Neutron Spectra 


Calculations of neutron-energy spectra based on measurements 
reported in Ringle's nee were repeated using our least-squares 
analysis. 

Ringle's experiment is briefly described as follows. Thresh- 
old detectors were exposed to neutrons from the Be *(d, n) Be reaction 
using 25- and 33-MeV deuterons from the 88-in, cyclotron at the 
Lawrence Radiation Laboratory. Threshold detectors were placedina 
ring 28 in. from the Be target and approximately 17 deg from the tar- 
get axis in the forward direction. Average deuteron-beam strengths 
for the bombardments were 2.37 pA for 25 MeV and 2.0 pA for 33 MeV. 


We used activities resulting from Ni (a, p)Co°’, Fe: ); 
g p 


56 48 a4 24 oT tn, evn Co? "(n, 2n) Co 


Mn°°, Ti*®(n, p)Sc*®, Mg@“(n, p)Na”*, Al pe 
reactions in recalculating neutron-energy spectra. 


and Nio’ (ny, anni 
Results of these calculations and previous measurements are shown in 
Figs. 18 and 19, Figure 18 demonstrates satisfactory agreement be- 
tween our method and previous measurement of similar neutron spectra 


by others, 


B. Neutron Spectra in a Thick Concrete Shield 
Exposed to a 6.2-BeV Proton Beam 


Our technique was next applied to the neutron spectrum inside a | 


thick concrete shield exposed to 6,2-BeV protons at the Lawrence Radi- &% 


ation Laboratory Bevatron. This experiment is reported elsewhere in 
detail. - Figure 20 shows a plan view of the shield, incident beam, 
and detector positions at which the neutron-energy spectrum has been 
solved, This spectrum measurement was complicated by two problems 
we had not encountered previously, One was that the variety of high- | 
energy secondary particles in the shield produced many radioisotopes 


in the detectors compared to the few isotopes resulting from Be "(d, n) Bi 
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Fig, 48, Measured neutron-energy spectrum from 
beam, §f 25-MeV deuterons on Be (17 deg from forward), 


Curve A has calculated by the step-function 


s been 
method, curve B by the polygonal tnethod, curve 

roblem C from emulsion data of Tochilin®‘ for 24-MeV | 

high- 1 deuterons on Be (forward), and curve D from | 
threshold-detector data of Heertje and Aten® for 

otopes 26-MeV deuterons on Be (forward), The curves 

7d: n)E of Tochilin and of Heertje and Aten were normal- | 


ized so that the area under the curves matches that 
under curve A. a 
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Fig. 19. Measured neutron energy spectrum from 33-MeV 
deuterons on Be (17 deg from forward), The two 
polygonal-method curves were calculated using iden- 
tical data but with different energy end points, 
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neutron sources, Gamma-ray spectra from detectors placed in this 
concrete shield were very complex, and some desired reactions, wi?’ 
(n, 2n) wi?! for example, were obscured by radioisotopes of similar 
y-ray energy. The second and more serious problem was that a signi- 
ficant fraction of the measured activity was caused by neutrons of 
energies in excess of 30 MeV, the upper limit of our spectrum solution 


technique. 


1. Correction for High-Energy Activation 

The fraction of activation caused by neutrons and other particle 
of energies greater than 30 MeV cannot be directly measured, Unlike 
the effect of cadmium on low-energy neutrons, there is no material 
with a large cross section which rises sharply below 30 MeV. An indi- 
rect method was used in which the shape of the neutron spectrum was 
estimated, and a high-energy threshold reaction used to determine the 


level or amplitude of that flux shape, 


a, Flux-Shape Estimate, Based in part, on the Monte Carlo calcu- 
lation of Metropolis et al., ae Moyer and Wallace" have calculated 
neutron-energy spectra for 450, 600, and 850-MeV protons incident or 
an extended thick object. From their number-density spectrum, n(E), 
a flux spectrum, v.n(E) = o(E), has been calculated (see Fig, 241). For 
convenience in calculations, the flux spectrum is approximated by an 
gt/2 shape above 30 MeV, For 30 < E < 1000 MeV the flux spectrum 


is assumed to be 


o(E) = 6(30)V 2, 


and for E> 1000 Mev, 


where $(30) is the particle flux at 30 MeV incm~“-sec”+-Mev~!. 
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Fig. 24, Calculated neutron-energy spectra for 850-MeV 
protons incident on an extended thick object, The 
solid line is an approximation to 6 (E). 
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Lacking appropriate calculations for 6.2-BeV incident protons, we have 
assumed that the neutron-energy spectrum has the same shape as for 


incident 850-MeV protons. 


b. Amplitude of High-Energy Flux, With rare exceptions, *+~45 


neutron-reaction cross sections have not been measured at high ener - 
gies, although a number of proton-reaction cross sections are fairly 
well known. In general, we have assumed that at the high energies, 
cross sections for formation of a particular radionuclide are the same 
for neutrons and protons. 

In polyethylene detectors placed in the concrete shield, we ob- 


and Cr tes 3p3n) Be. The cross sec- 


served Bee from Na een a2n} Be 
45,46 


i 


tion for oats 3p3n) Be! is a well-determined function of energy 


and is shown in Fig. 22. The threshold for this reactiecn is 38.5 MeV... 


Under our assumptions, the Be! activity 1s expressed as 


| os 4000 
_ _ : = {30 
A, -/ $(E) o,(E)dE = (30) f == 9,(E) dk, (417) 
: 0 0 
i and therefore, 
: A,' 
$30) 7000 


Having calculated (30), we use this same procedure to calcu- 
late the activation of other detectors caused by neutrons of energies 


above 30 MeV, 





c, High-Energy Cross Sections, Figures 23 through 27 show neu- 
tron cross sections below 30 MeV and proton cross sections above 30 
MeV. Sources of experimental proton cross sections are numerous, 
and are not identified in the following figures, The smooth curves were 
drawn to define o(E) used in the calculation of 
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Fig, 23, Excitation function for Mg" 24 P p)Na“4 . Data 
points are shown for Mg(p, -)Na@ See text for an 
estimate of the high-energy cross section, 
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Fig, 24, Excitation function for Al*"(n, a)Na“*, 
points are shown for Al“? (p, 3pn)Na** 
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Fig, 25, Excitation function for Fe tn, p)Mn°°, A data 
point is shown for Fe(p, -)Mn° - See text for an 
estimate of the high-energy cross section, 
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Fig. 26, Excitation function for Co? (n, 2n Gao. Data 


points are shown for Co”?7(p, pn)Co°°, 
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A literature search did not reveal any experimental cross- 
section values for ors formation by neutrons or protons above 30 MeV 


on nickel, The high-energy cross section was estimated as follows: 


58 


At 396 MeV, the cross section for Ni(p, pxn) Gar? is 45mb. Naturally 


occurring nickel has five isotopes: ni?3 67.8% abundance; nie? 


A 
26.2%: ni?! 12%; wie’, 3.7%; and nio4 


assume that Ni? tp, 2pn) Co?” is the predominant reaction, and that 


, 1.2%. At high energy, we 


others resulting in Cae [| for example, NiO %p, 2p3n) Co”7 may be 
neglected, From this assumption it follows that the Nico. 2pn) Ger? 
cross sectionis 45X1000/0.678 = 66 mb. For formation of Go? the 
predominant reaction at high energies is assumed to be Ni tp 2peyGe’. 
and other reactions resulting in Ger are assumed to be negligible. 
Next we assume that the cross section for Mio Gp: 2pn)Co?® is also 66mb, 
It follows that the formation cross section for €or" from natural Ni 
(per NiO” nucleus) is 66X0.262/0. 678 = 25 mb. As before, we assume 
that the high-energy reaction cross section is the same for neutrons 
and protons, 

Experimental high-energy cross-section data is available for 
Mg(p, saynal’. ayine 
an average value of about mb above 100 MeV. Isotopic abundances in 
Mg are Mg““, 78.7%; Mg~°, 10.14%; and Mg@©, 14.3%. Per Mg@* 


nucleus, this cross sectionis 6/0.787 = 7.6. Protons may produce 


A plot in Fig. 23 shows that this cross section has 


Na“ in Mg in the following reactions: 


25* Natt e p 


1. Mg@*(p, p')Mg 
23 Mg? “(p, p'yMg7°* > Na“4 +ptn 
3, Me“, a) Ale’ a as? a 2p: 
Neutrons may produce Na“4 in Mg in the following reactions: 


4, Mare Ga: Nae 


5. Mg~”(n, aioe?” > Na + p 

6. Mar x: Siac” > Nace +n 

te Me: n! yMg*o* > Na“4 +ptn 
8. Mg~"tn, p)Na*°* > Na“4 + 2n 














Based on measurement of the cross section for roger pipe at 5.4m. 
at 370 MeV, a reaction 4 is estimated to have a cross section of 5 mp. 
Reactions 5 plus 7 plus 8 are assumed to have the same cross section 
as 1 plus 2 plus 3. To account for reaction 6 for which there is no 
similar proton-induced reaction, the total for numbers 5 through 8 is 
assumed to be 1.5X7.6 = 14.4 ms. It follows from these assumptions 
that the high-energy-neutron reaction cross section per Meg** nucleus 
is 14.4 +5 46mh. 

There are no experimental-cross-section measurements for 


r 


Mar’ production by high-energy neutrons oniron. The 340-MeV 
cross section is 0.59 mi. for Fe(p, speci” and 0,24 mi for Fe(p, xn) 
Ge?” 28 Naturally occurring iron has four isotopes: Heoe, 5.8% 
abundant; Fe°®, 91.7%; fee 2.2%; and peo 0.3%. High-energy 
protons on iron may. produce Mn°° by the reaction: 


st 5 
73 56 i 


1, Fe (p, p')Fe?’ + Mn’? + p. 


Here Fe?? is neglected, being in very small quantity. High-energy 
neutrons may produce Mn” by: 


2. Fe wae p)Mn°° 
6 


3. Bae (a; n')Fe!* > Mn? + p 
4, Fe? (n, p)Mn°” ° 


ste 
os 


Sian +n, 


The cross section for reaction 2 is estimated as 0.24 mb from the 


2 56 fs 
cross section measured for Fe(p,xn)Co” -. The cross section of re- 


actions 3 plus 4 1s estimated to be twice that of 1, From these assump- 
56 


tions, it follows that the high-energy-neutron cross section for Mn 


ie 
production, per we?? nucleus, is 2X0.59 + 0.24 = 1.4mb, 


d. Calculation of High-Energy Activation, For each location in the 
shield where the Ge te a2niBe” reaction was observed, the calculated 


activation due to particles of energy greater than 30 MeV is 


1000 


Peasy =f $(E)o,(E)dE = $(30) 4 3 o(E)dE, (18) 
* 30 30 
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From ‘Ago, = A, & Amax, » where A, is measured and 


ae is calculated above, we can solve the basic set of equations for 


6(E) in the previous manner. Because the high-energy flux shape and 


> 


cross sections are both estimates, the calculated activation, Bohs 


is subject to a large error, This error is taken to be of the order of 


A ‘ 
max 
A’; in one case, it was 23.5%. 


For most of- our. reactions Ba was between 1 and 10% of 


2, Measured Spectra 

Figures 28 through 32 show neutron spectra determined at a 
number of locations in the shield, Note that these spectra decrease 
with increasing energy, as expected. The high flux values below 10 
MeV are consistent with emulsion measurements of Bevatron stray- 
neutron spectra, 27 The spacing of reaction thresholds is such that a 
peak in the few MeV region could not be resolved in these calculations, 

At each location except 4' - 1', five reactions were used: 


Nie ern, p)Go"®; ee OG, p)Mn-°; A ajNacss Mg“ *(n, p)NaZ*; and 


Co°%n, 2n)Co>8, 127 126 


At location 4' - 1', the reaction I {n, 2n)I was 


added to these five. 


3, Evaluation 

Note that the spectral solutions are stable, in that they are not 
oscillatory. This was assured by limiting the number of variables to 
three, for example, when five reactions were used, Such a problem 
would be over-determined if solution were attempted by matrix inver- 
sion; however our least-squares technique makes an optimum selec- 
tion of three flux variables to fit five activities. 

If experimental activity measurements are subject to large er- 
ror, as was the case in our experiment at the Bevatron, oscillatory 
solutions can result from our least-squares technique. Solutions can 
be made stable by limiting the number of variables, Figure 33 dem- 
onstrates this behavior for spectral solutions at location 8'- 4! in 


the Bevatron experimental shield, It will be noted that the step-func- 


tion solution shows more stable behavior than the polygonal solution, 





Comparison between measured activities and activities calculated 
from the spectrum solution indicate that measured activities contained 
errors as large as 25% for threshold detectors exposed at location 

8! - 4', 





-53- | 


pie) (n em72 sec! Mev~') x 10° 





MU-35006 


Fig, 28, Measured neutron-energy spectrum in the 
Bevatron shield at location 4'-4', using step func- 
tion and polygonal representations, One polygonal 
curve (long dashes) is based on measurements of 
five reactions, The other polygonal curve (short 
dashes) is based on measurements of six reactions, 
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Fig.29, Measured neutron-energy spectrum in the Bevatron 
shield at location 4'-3', using step-function and poly - 
gonal representations, 
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Fig. 30. Measured neutron-energy spectrum in the Bevatron 
shield at location 8'-2', using step-function and poly- 


gonal representations, 
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Fig, 34. Measured neutron-energy spectrum in the Bevatron 
shield at location 8'-4', using step-function and poly- 


gonal representations, 
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Fig, 32, Measured neutron-energy spectrum in the Bevatron 
shield at location 12'-3', using step-function and poly- 
gonal representations, 
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Fig, 33. Measured neutron-energy spectrum in the Bevatron 


shield at location 8'-4', The same five reactions were 
used for all spectra, Better stability is observed for the 
step function, The two polygonal curves in the bottom 
graph show that an unstable solution has no particular 


pattern, 
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VI. SUMMARY 


A method of measuring neutron spectra in the 2,5-to-30-MeV 
range, using threshold detectors has been improved and tested in 
application. Threshold reactions (n, p), (n,a), and {n, 2n) are used 
where the residual nucleus emits a y ray and where neutron cross 
sections have been measured as a function of energy. Detectors are 
counted with a Nal(T1) scintillation crystal and resulting y spectra 
are analyzed by a computer program (SUPER-3). Two new computer 
programs (FLUXSPOS and FLUXPPOS) use least-squares analysis 
to solve the neutron spectrum based on input ‘of activation and re- 
action-cross-section data, 

The y-spectrum-analysis method is a previous development 
by Dr. John Ringle* and is briefly described earlier in this report. 

Our improved method of neutron-spectrum analysis uses ex- 
perimentally determined cross sections in place of cross-sections cal- 
culated by continuum-model theory. Data for cross-section measure- 
ments are taken from current literature. | 

In measuring the spectrum of a very-high-energy neutron 
source, the fraction of total activity caused by neutrons of energies 
below 30 MeV is estimated. 

Detector activity is related to cross section and flux through 
Eq. (5). Once a form of the flux spectrum has been chosen (for ex- 
ample step function or polygon), this integral equation set may be 
solved by matrix-inversion techniques. However, such techniques 
result in oscillatory solutions where activity measurements contain 
substantial experimental error. A new method of solution is usedin 
which least-squares analysis chooses flux variables such that an op- 
timum fit is made between measured activities and activities calcu- 
lated from the cross sections and flux solution, Using our technique 
we were able to measure spectra satisfactorily where experimental 


errors were unavoidably large. 























The neutron spectrum of Be (d, aye? was measured using 25- 


MeV deuterons and compared with other measurements of that spec- 
trum using 24-MeV and 26-MeV deuterons, Comparison indicated 
that our technique was satisfactory. The spectrum was also measured 
for this reaction using 33-MeV deuterons, 

Our technique was applied to the neutron spectra at several 
locations in a thick concrete shield exposed to 6,2-BeV protons, In 
this case the activation caused by neutrons of energies greater than 


30 MeV was estimated by: 


14. estimating the high-energy flux shape; 

2, calculating the magnitude of that flux by measuring Be! 
activity in carbon (38. 5-MeV threshold); 

3, estimating high-energy-neutron cross sections on the basis 


of measuzved proton cross sections. 


Although the complex y-ray spectra in these detectors limited accuracy 
in determining activities, we were able to measure neutron spectra 


successfully with our new least-squares technique, 
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VIl. CONCLUSIONS 


Our y-spectrum analysis method, SUPER-3 developed by 
Ringie* is quite satisfactory, It is more suitable than other methods 
for this application because of the accuracy with which it analyzes 
y-ray spectra containing a number of unknown radioisotopes in addi- 
tion to the one being analyzed. A very time consuming determination 
of the identity and response functions of these unknown radioisotopes 
would be necessary if least-squares y-spectrum analysis were used, 

New measurements of cross section for (n, p), (n,a), and (n,2n) 
threshold reactions have been published and show sufficient resolution 
that excitation functions can be determined by drawing smooth curves 
through the data. Continuum-model calculations of cross sections do 
not dependably agree with measurements, 

" Two computer programs, FLUXSPOS and FLUXPPOS, suc- 
cessfully compute neutron spectra from activity measurements and 
cross sections, These programs are usable where activity data con- 
tain large experimental errors. Toa great degree, these two pro- 
grams avoid unstable and misleading solutions which other methods 
tend to give. | 

This improved technique is applicable to neutron spectra in the 
Bevatron shield, when suitable correction is made for activation by 
high energy particles (>30 MeV). It should be applicable to stray- 
neutron spectra around the Bevatron. 

The resolution of this method should be improved by using more 
reactions than were used in spectra reported here (five to seven), It 
should be of interest to investigate the resolution of this method in 


solving for a known, artificial spectrum with known inaccuracies in 


activities, 
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APPENDIXES . 


A. Equations for Detector Activity 


Our computer program, SUPER-3, determines the photopeak 
area, ae by fitting a Gaussian distribution to the corrected gamma- 
spectrum data (Refer to paragraph II. B.1.) The activated threshold 
detector is analyzed for y radiation and is referred to below as the 


source, This program computes the source disintegration rate by 


solving 
as p 
An e Bpa(1-BT)’ v20) 
Pp 
where 
Ap = source disintegration rate in disintegrations per second, 
NS = photopeak area in counts in peak per second, 
Ep = peak efficiency in counts in peak per gamma emitted by the 
source, 
B = gamma branching ratio in gamma emitted by the source 
per disintegration, 
a = correction factor for source self-absorption, and 
BT = block time = fraction of the time the PHA does not accept 


pulses, 


The correction factor, a, for source self-absorption may be defined 


further as 


4 Pa: 
secs, (20) 
where 
T = linear y attenuation coefficient in em, 


t = source thickness in cm, 











The source disintegration rate, An may be corrected to 


saturation activity per target nucleus, A', by using 


Ag At 
i 21 
: as a oi 24) 
N(i-e je d-e 
where 
: N = total number of target nuclei, 
» = decay constant for the residual nucleus in min , 
| 
| ty = irradiation time in minutes, ‘ 
| t., = time between the end of irradiation and the beginning of 
y counting in minutes, 
to = time for y counting in minutes. 
If the detector is in the form of a pure element, we have 
- w No r w Na # 
ep oars © age 
Ww p 
where 
w = detector weight in grams, 
Ny = Avogadro's number, 6.02 «1073, 
A, = atomic weight of target element, 
r = abundance of the target isotope in the target element, 
Ny, = atomic density in cem™?, 


: : ~3 
p = mass density in grams -~cm -, 


B. Reduced Experimental Data 


The following table shows A390, [Eq. (5)],the detector-saturation activity caused by neutrons of less 
than 30-MeV energy. These saturation activities are computed from Eq. (21) and include correction 
for high-energy activation calculated by Eq.(18). The units are disintegrations per second per target 
nucleus times 10°~. 


Reaction Experiment 
Deuterons on Be Bevatron shield 

P.O > “s 

(d' energy indicated) (Locationindicated) 

25 MeV 33 MeV. Ala 4\ -3! gt -2' gi 4! 421-3! 
Ni? ®(n, p)co? 754,470 935,770 769,170 67,818 72,042 414,069 5,369.8 
Be (a p\Mn-” 70,845 55,184 5,068.5 4,868.0 804.17 334.98 
Ti’ (appise” 47,388 56,514 
Me“ ta, piNa’” 454,980 254,270 86,144 8,313.5 8,431.7 4,970.9 719.64 
Al“'(n,ajNal? 88,259 140,020 55,864 4,900.3 5,234.6 1.188.6 436.44 ci 
fee Sao 575,150 
Co (a, 2njGor” 364,930 776,860 357,970 32,527 34,706 7,990.4 3,145.0 
Mi?’ ann?! 24,338 54,438 
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